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Abstract 

A novel robust mechanism for the generation of "trapping states" is shown to exist in the 
coupling of a two-level system with an oscillator, which is based on nonlinearities in the 
laser-induced vibronic coupling. This mechanism is exemplified with an ion confined in the 
potential well of a trap, where the nonlinearities are due to Franck-Condon type overlap 
integrals of the laser waves with the ionic centre-of-mass wavefunction. In contrast to the 
coherent trapping mechanism known from micro-maser theory, this mechanism works also 
in an incoherent regime operated by noisy lasers and is therefore much more robust against 
external decoherence effects. These features favour the incoherent regime, in particular for 
the preparation of highly excited trapping states. 



1 Introduction 

Micro-maser theory predicts a coherent regime where so-called trapping states are 
created in the microwave cavity field, by injection of electronically prepared Ry- 
dberg atoms with well-defined velocities [1]. They are stable quantum states of 
the cavity field being unchanged by further interactions with injected atoms. The 
trapping states are photon-number states that represent a discrete, fixed number of 
cavity photons [2] . The preparation of trapping states with a pre-described photon 
number is performed by adjusting the atom-field interaction time, i.e., the velocity 
of the injected Rydberg atoms, to certain values. This mechanism crucially depends 
on a narrow velocity distribution of the atoms and is also very sensitive to any deco- 
herence effects, such as cavity decay and thermalization via background radiation, 
or thermal motion of the cavity mirrors. Another limitation is the occurrence of 
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(rare) two-atom events, that destroy the trapping states. This sensitivity is due to 
the fact, that the trapping mechanism relies on the coherent features of the Jaynes- 
Cummings dynamics [3]. Therefore, initially not a perfect photon-number state 
but sub-Poissonian photon statistics have been realised [4]. Later on, signatures of 
trapping states with small photon number could be successfully demonstrated [5]. 
When preparing trapping states of relatively large photon numbers, however, the 
above mentioned detrimental influence of decoherence poses an essential problem 
for the application of the coherent trapping mechanism. 

Alternatively, a coherent Jaynes-Cummings type trapping mechanism can also be 
realized in the motion of a single trapped ion [6,7,8,9]. Here the harmonic vibration 
of the trapped ion plays the role of the micro-maser cavity field, and the interaction 
of the injected Rydberg atoms with the micro-maser field is implemented by a 
combination of optical pumping and optical excitation of a vibrational sideband. 
As in the case of the micro-maser, this trapping mechanism for the vibration of a 
trapped ion also depends crucially on a coherent Jaynes-Cummings dynamics [7]. 
For a realistic experimental situation with phase-fluctuating laser fields obviously 
the coherent trapping mechanism is disturbed. In that case it has been shown that a 
sub-Poissonian, binomial statistics emerges for the vibrational quantum number of 
the trapped ion, that reveals a noise level being half of the classical shot noise [9]. 

The crucial point for realizing Jaynes-Cummings type trapping states in the motion 
of a trapped ion seems to be the need of a coherent time evolution without any 
decoherence effects. If one could use incoherent dynamics to prepare the desired 
states, all the above-mentioned problems could be avoided. Thus the question arises 
how to realize an incoherent dynamics that also gives rise to appropriate trapping- 
state conditions. 

In this paper we propose a novel mechanism for the generation of "trapping states" 
in the vibration of a trapped ion that is based on nonlinearities occurring in the 
laser-induced vibronic interaction [10]. Dephasing effects emerging for example 
from phase-fluctuating laser fields or spontaneous emissions do not disturb this 
mechanism. Therefore, this incoherent regime is realistic and is also experimentally 
more feasible than the coherent dynamics needed for micro-maser-type trapping 
states. Moreover, the incoherent trapping mechanism may be less sensitive to other 
types of externally induced decoherence. 

In Sec. 2 the laser-excitation scheme and the corresponding interaction Hamilto- 
nians are discussed. In Sec. 3 we briefly summarise the realization of a coherent 
micro-maser-type dynamics, and in Sec. 4 the incoherent regime is introduced. A 
trapping mechanism that works also in the incoherent regime is then discussed in 
Sec. 5. Finally, a summary and some conclusions are given in Sec. 6. 
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Fig. 1. Laser-excitation scheme for the preparation of trapping state. During the opti- 
cal-pumping interaction (a) a laser (Rabi frequency Q p ) together with spontaneous emis- 
sions (rates 712) prepares the system in state |2). In the sideband interaction (b) the first 
vibrational sideband (vibration frequency v) is driven (Rabi frequency fi s ) either directly 
(solid arrow) or in a Raman configuration (dashed arrows). 

2 Laser excitation scheme for the trapped ion 



The centre-of-mass coordinate of a single ion bound in a radio-frequency trap, such 
as for example a Paul trap, is subject to a dynamics that can be described to good 
approximation as a 3D harmonic oscillation with three characteristic vibrational 
frequencies along the principal axes of the trapping potential. The geometry of the 
propagation of externally applied laser fields can be configured in such a way as to 
influence only the ion's centre-of-mass motion in the direction of one principal axis, 
that is associated with the trap frequency v. Thus, one may consider the system as 
being one-dimensional, due to the natural decoupling from the remaining two other 
degrees of freedom [11]. 

In what follows, we always consider the case of quasi-monochromatic, collimated 
laser beams, whose wave-vectors are parallel to the chosen principal axis of the 
trap, that we may specify by the centre-of-mass coordinate x. In the case where 
laser-atom interactions are implemented by two-photon induced Raman transitions, 
the wave-vector of the beat note of the required two laser beams shall point along 
this axis. 

The laser-excitation scheme shall consist of cycles, each being implemented by 
two interactions, an optical-pumping interaction, and a sideband interaction. The 
optical pumping is intended to invert the populations of the two states |1) and |2), 
cf. Fig. la. It is implemented via the laser-driven electric-dipole transition 1 1) |3), 
and the spontaneous decay on the electric-dipole transition |3) — > |2). The equation 
of motion, governing the dynamics of the vibronic density operator of the ion, g, is 
then given by 
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-^[H +H p (t),g]+^ 1 (AuK^g] A 31 - \A 33 g-\gA 3 ^ 

+ 72 (A 23 n 2 [Q]A 3 2 - \A 33 g- ±qA 33 ) . (1) 



The operators A\j = \ induce transitions between the ion's electronic states \j) 
and |z) and 71 and 72 are the spontaneous decay rates of state |3) decaying to states 
|1) and |2), respectively. The free Hamiltonian H reads 

3 

H = ^ Q^a + ^2 huiAu, (2) 
i=i 

where a and a) are the annihilation and creation operators of vibrational quanta, 
respectively. 

The laser interaction driving the optical pump process reads 

if p (t) = i^ p i 3 ie-^ i + H.a., (3) 

where Vt p is the Rabi frequency denoting the coupling strength of the pump laser 
beam with the electric-dipole moment of the ion. The pump-laser frequency uj p 
shall be near resonant with the electronic transition frequency to 3 i of the states 
|1) and 1 3). For this interaction the laser beam is supposed to propagate along a 
direction perpendicular to the motional direction x under consideration, to avoid 
coupling of the electronic and vibrational degrees of freedom in x direction. 

The photonic recoil acting on the centre-of-mass motion of the ion during the spon- 
taneous emissions is given by the operator TZi (i — 1, 2) 

1 

Ki[Q\= ( dsw(s)e lsr] ^ +h) Qe- isr >^ +h \ (4) 



with w(s) = |(l+s 2 ) describing the dipole radiation characteristics and the Lamb- 
Dicke parameters r]i and 7] 2 correspond to the spontaneous transitions |3) — > \ 1) and 
1 3) — > 1 2), respectively. For simplicity we neglect, for the moment, the vibrational 
scattering during the spontaneous photon emissions, i.e. we choose TZi[g] — > g, 
though later these effects will be taken into account to full extent numerically. 

The stationary solution of the optical pumping process is approximately reached 
after a sufficiently long time of laser interaction. In terms of electronic-state matrix 
elements = (i\g\j) of the density operator, that are still operators acting on the 
vibrational degree of freedom, the stationary solution reads 

Pi2(r p ) = p 23 (T p ) = /3 3 i(r p ) w (5) 



4 



PnOp) = /5 3 3(Vp) « (6) 
p 22 (r p ) = U (r P ) [Pii(0) +p 22 (0)] Ul(r p ). (7) 

Here we assume that initially, due to relaxation, the population in state |3) is zero, 
i.e., P33(0) = 0. Thus, all population from state |1) merges with the population of 
state 1 2), while the ion is subject to oscillation in the trap potential. 

The subsequent sideband interaction is implemented either by direct or by Raman 
excitation, resonant to the red vibrational sideband of the weak transition 1 1) <->■ |2), 
cf. Fig. lb. Its interaction Hamiltonian reads 

H s {t) = \mA 21 e 1 ^-^ + H.a., (8) 



where k and to are the projection of the wave- vector on the x axis and the frequency, 
respectively, either of the laser field for direct excitation, or of the beat note of the 
Raman lasers. Moreover, the Rabi frequency f2 denotes the coupling strength of this 
transition. For the case of resonance with the red vibrational sideband, uj = uj 2 i — v, 
this laser interaction induces resonant vibronic transitions between states |2)|n) 
and |l)|n+l), where |n) are energy eigenstates of the centre-of-mass vibration 
in the trap potential. For this case the Hamiltonian (8) simplifies, in rotating-wave 
approximation and in the interaction picture, to a nonlinear Jaynes-Cummings type 
Hamiltonian [10], 

H s = ^hQ s A 21 f(h; r))a + H.a. , (9) 

where f2 s = ir)Q is the Rabi frequency on the red vibrational sideband in the 
Lamb-Dicke limit. The effects that arise in the regime beyond the Lamb-Dicke 
approximation, i.e., with rather large Lamb-Dicke parameter rj are described by 
the excitation-dependent operator function f(h;rj) with h being the number op- 
erator of vibrational quanta in x direction. Note, that for a Raman configuration, 
changing the laser-beam propagation geometry varies the projection k of the beat- 
note wave-vector on the x axis. Thus, the Lamb-Dicke parameter can be tuned up 
to rather large values. 

The operator function f(h; 77) depends solely on the number of vibrational quanta 
and can be given by a normally ordered expression as 



00 2n 

/( ft;l) ) = e-^|; ( _ ir _l_ a »»a- (io) 

= : (^VfMy 1 J x (2r/v / ata) : e^ 2/2 , 

where : : denotes normal ordering, and Jk(z) are the Bessel functions of integer 
order. Clearly these functions introduce a dependence of the laser-ion coupling on 
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the vibrational excitation of the ion. They can be understood as analogy to Franck- 
Condon overlap integrals, however, in momentum space due to the recoil of ab- 
sorbed and emitted photons. 

For this approximation, i.e. neglecting spontaneous recoil effects, in each cycle of 
optical pumping and sideband interaction, the mean vibrational excitation will be 
increased or at least held constant. The dominant processes in one cycle can thus 
be given by the cascade 

pump v sideband \^)\ n ) 



\2)\n) 

\2)\n) |l)|n+l> 



thus vibrational transitions |n) — ► |n) and \n) — > |n+l) are performed in each cycle, 
leading to a net increase of the vibrational excitation. 



3 Coherent time evolution 



In the absence of decoherence effects the time evolution for the sideband interac- 
tion is described by a unitary time-evolution operator U (t). Thus, starting with the 
vibronic density operator g of the trapped ion [cf. Eq. (7)] 



(t+T p ) = A 22 U (r p ) [g u (i) + g 22 (t)} Ut(r p ), (11) 



at the time t+r p , right after the optical pumping, the density operator reads after an 
interaction time t s of the sideband interaction, 

g(t+T p + r s ) = U s (t s ) g(t+T p ) ^(r s ), (12) 



where U s (t) = U (t) exp(—iH s t/h) and U (t)=exp(—iH t/h). 

Starting from an initial density operator at time t k we obtain the density operator at 
the time t k+ i —t k +r p +r s , after a complete cycle of pump and sideband interaction, 
by insertion of Eq. (11) into (12) as 

g(t k+1 ) = U s (t s ) |2> U (r p ) [(l\g(t k )\l) + (2\g(t k )\2)] t)J(r p ) (2\ fr s (r s ) (13) 



The vibrational number statistics is obtained from the density operator as 

P n (t) =Tr [g(t)\n)(n\], (14) 
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Using Eqs (13) and (14) straightforward calculation [13] results in a recurrence 
relation for the probabilities of vibrational quantum numbers 

P n (tk+i) = w coh (n\n-l) P n -i(t k ) + w coh (n\n) P n (t k ), (15) 



where the probabilities for excitation or survival are characterised by a coherent 
sideband-interaction and read 



^coh(n+l|n) = sin' 
w co h(n\n) = cos 



f(n;rj)y/n+l 



h\to B \T a f(n;ri)y/n+l 



(16) 
(17) 



where we have used the function f(n;rj) = (n\f(h; rf)\n), which can be derived 
from Eq. (10) as Laguerre polynomial 



/ (71577) = (ra+1)' 



(18) 



Equation (15) describes the change of the vibrational statistics depending on the 
number of cycles performed. This recurrence relation contains the familiar trapping 
mechanism known from micro-maser theory [1]: The trapping state |n ) is reached 
when the excitation rate (16) from state \n ) to \n + 1) vanishes, so that all the 
population accumulates in state \n ). The condition w; coh (ri + l|?i, ) =0 thus leads 
to the condition for \n ) being a trapping state: 

\fl s \r s f (no; rj) \/Vio + l = 2-nm, m G Z. (19) 



In the Lamb-Dicke limit the nonlinear coupling function in Eq. (19) reduces to 
lim^o f{n; rj) = 1, and Eq. (19) reduces to the familiar trapping-state condition for 
the micro-maser case [9]. This case is characterised by a non- vanishing function 
/(no; rj), upon which the condition for the interaction time results 

2irm ^ ^ 

\n s \f(no;r])^/no + l , 



This type of trapping states is generated by a coherent mechanism, that consists of 
complete Rabi cycles at the trapping- state condition, leading to a locked vibrational 
quantum number. 

Eq. (19) on the other hand, contains also a quite distinct type of trapping condi- 
tion, namely that one where the Franck-Condon-type nonlinearity in the laser-ion 
coupling strength leads to a vanishing coupling, i.e., if 

/(n ;77) = 0. (21) 
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Fig. 2. Coherent transition rates w co \ 1 (?i+l\n) for (a) ?] = 0.1, |il s |r s = 1.149 [trapping-state 
condition (20) for no = 50], and for (b) 77 = 0.268 [trapping-state condition (21) for uq = 50] 
and |r2 s |r s = 7r/2. 

This case does not correspond to complete Rabi cycles, but represents the break- 
down of the laser-ion coupling mechanism, due to non-overlapping trap eigenstates 
in momentum representation, that are shifted off each other by the differing amount 
of the effective photonic momentum hk. This novel type of trapping states may be 
set up in both a coherent or incoherent way. 

Let us illustrate the qualitative differences between the two types of trapping states 
for a coherent sideband interaction and consider the two different ways of generat- 
ing a trapping state \n ) with no = 50 vibrational quanta. The well known way is to 
use the standard trapping- state condition (20) to fix the pulse area \Cl s \r s of the laser 
resonant with the vibrational sideband for a given Lamb-Dicke parameter n. The 
alternative way of generating this trapping state in a coherent fashion is given by 
adjusting the laser-beam geometry for sideband interaction (in a Raman configura- 
tion) in such a way as to tune to a Lamb-Dicke parameter that fulfils condition (21). 
The laser-pulse area is then chosen in such a way, that lower-lying trapping states 
due to complete Rabi cycles are avoided. 

The improvement is illustrated when comparing the coherent transition rates w co h(n+ 
l\n) in the two cases, cf. Fig. 2. The dip of vanishing transition rate around the trap- 
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ping state \n = 50) is sharper for the rate shown in part (b), which is due to the 
nonlinear modification encoded in the function f(n; 77) when 77 is not too small. As 
a consequence, when the quantum state approaches the trapping state \n = 50), the 
dynamics is significantly slowed down in case (a) since the rates u; coh (n + l|n) are 
smaller (for n = n — 1, n — 2, n — 3) as in case (b). In both cases the transition 
rates may be suppressed at certain vibrational numbers below the trapping state. 
These coherent effects emerge when the ion performs nearly complete Rabi cycles 
on the electronic transition at certain values of the integer m in the trapping- state 
condition (20). Such subharmonic resonances degrade the generation of the trap- 
ping state, since they substantially decrease the flow of population into the trapping 
state. 

The sensitivity to subharmonic resonances seems a handicap to both coherent ways 
of achieving a trapped state, in addition to the fundamental problem of preserving 
coherence. The novel scheme, however, evades both problems, since it is applicable 
in an incoherent way. 



4 Incoherent dynamics 

The trapping mechanisms described above, being based on coherent temporal evo- 
lution, are disturbed by any decoherence, as for example the phase fluctuations of 
the Raman lasers, or spontaneous emissions. It would be favourable to have a trap- 
ping mechanism at hand, that does not depend on the coherence of the underlying 
physical processes and is therefore not disturbed by such effects. Then, a trapping 
state might be reached even under the influence of various decohering perturbations 
and would be more robust than the micro-maser-type trapping states. 

We consider the sideband interaction under the influence of decoherence, i.e. elec- 
tronic dephasing, such as produced by phase-fluctuating lasers, spontaneous emis- 
sions, or collisions with background vapour. The total amount of phase fluctuations 
will be described by a dephasing rate 7 that gradually destroys the coherence be- 
tween electronic levels |1) and |2). For simplicity we do not consider electronic 
dephasing mechanisms that depend on the vibrational excitation. Such effects may 
appear due to laser fluctuations [14], spontaneous Raman processes [15,16], or 
combinations of both [17]. Moreover, for the analytic treatment we neglect the 
spontaneous recoil effects. These however will be included in the complete numer- 
ical solution shown later. 

The master equation for the sideband-coupling interaction, including the decoher- 
ent phase fluctuations reads in the interaction picture as 

h = - l j\H s ro\-\ (in Q A22 + A 22 q i n ) , (22) 
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where the second term accounts for the electronic dephasing [18,19]. The equations 
of motion for the electronic populations are derived from Eq. (22) as 



$11 = - 2 fa* s &Q 21 -n B Q 12 b), (23) 

&2 = -^(fi B 6ei2-^&i6 t ) > (24) 
where we have made use of the nonlinearly deformed annihilation operator 

b = f(cra; rj) a. (25) 

The time evolution of the electronic coherences, on the other hand, is modified 



02i = _ 7r£2i - 77^s (bgn - 822b) . (27) 



£12 = ~\ Q12 ~ (b ] h2 - Qii , (26) 

7 » i. 
2 621 ~ 2 ! 

Let us now consider the incoherent regime where the dephasing rate 7 is larger than 
the Rabi frequency on the vibrational sideband: 7> |f2 s |. Nevertheless, we consider 
the case of well-resolved sidebands, v > 7, and thus 7 is in the range given by 

|Q„| < 7 < v. (28) 

Then the rotating- wave approximation with respect to v for the derivation of Eq. (9) 
is valid. For typical vibrational frequencies of u/2ir ~ 1. . .10 MHz and Rabi fre- 
quencies of |fi|/27r ~ 10. . .100 kHz, the dephasing rate is supposed on the order 
of 7/27T w 100 kHz . . .1 MHz. In this incoherent regime adiabatic elimination of 
the electronic coherences can be performed by neglecting the time derivatives in 
Eqs. (26) and (27) and solving for the electronic coherences. Inserting the resulting 
electronic coherences into the equations of motion for the electronic populations 
(23) and (24) we obtain 



S11 = — (#822 b - \#b Q U - \hx Pb) , (29) 
'SeuSt-iSSt&a-i&zSSt). (30) 



022 
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These equations can be decoupled and solved by use of the dressed-state coeffi- 
cients 

C mn = (m\Qu\ n ) ± (m-l\g 2 2\n-l). (31) 
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From Eqs. (29) and (30) we obtain the decoupled differential equations for the 
coefficients C^J (t) as 



2 

C7(±) = _I 

^— ^ <rr> <n A 



mn 4 



^Jlm-MT^ln-M C£>, (32) 



where the excitation-dependent damping coefficients are given by 

7n (r / ) = 2 7S (n+l)[/(n;r / )] 2 , (33) 

and the saturation parameter reads s = 2|f2 s | 2 /7 2 . Using Eq. (18) the nonlinear 
damping coefficients can be given in terms of Laguerre polynomials as 

ln(v) = ^\LH\v 2 )\ 2 e-* 2 . (34) 

72 + 1 

Note, that the laser-ion coupling nonlinearity appears only in these damping coef- 
ficients, and that in the Lamb-Dicke approximation they become linear with n: 

lim 7n (r/) = 2 7 s(n+l). (35) 
n— »o 



The vibrational number statistics P n can be easily obtained from the matrix ele- 
ments C^l as follows: 



Thus only the diagonal matrix elements C^' are needed, which obey the equations 
of motion 

CS>=0, CtJ = -ln-i(v) CL ] . (37) 



Since the incoherent dynamics of the sideband interaction follows the optical pump- 
ing, the initial conditions for the density-matrix elements can be written as 



(n\g u (t k + T p )\n)=0, (38) 
(n\(hi(t k +T p )\n) =P n (t k + T p ) = P n (t k ), (39) 

where again t k denotes the time after a complete cycle, and r p is the pump time. 
Using these initial conditions one obtains, via Eq. (31), the initial conditions for the 
coefficients and the solution of Eq. (37) results as 
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CM(t k+1 ) = P n ^(t k ), (40) 
Ci~\t k+1 ) = -Pn-xih) expPy^fo) r B ], (41) 

Inserting Eqs. (40) and (41) into Eq. (36) one finally obtains a modified recurrence 
relation for the number statistics after subsequent cycles in the incoherent regime 

Pn(h+i) = w inc (n\n-l) Pn-\(tk) + w inc (n\n) P n (t k ), (42) 

with the transition rates in the incoherent regime of sideband interaction, 

w- mc (n+l\n) = \ {1 - exp[-7„(?7) r s ]} , (43) 
Wine ON = 1 1 1 + exp[-7 n (r/) r s ]} . (44) 



5 Incoherent trapping mechanism 

Let us now study in more detail the rates Wi nc (n+l\n), that determine the increase 
of the mean vibrational excitation. If those rates are non- vanishing, the populations 
of vibrational quantum numbers will be shifted to the next higher lying quantum 
numbers in each cycle, leading to an overall increase of the mean vibrational ex- 
citation. If, however, Wi nc (n + l\n ) is zero for some vibrational number n Q , the 
recurrence relation (42) has a cutoff at that number no. That is, all populations that 
were initially at vibrational quantum numbers below n will be accumulated after 
some time in the vibrational state |n ). Clearly, populations in levels higher than n 
will be shifted to still higher quantum numbers in each cycle. 

Thus two distinct dynamical regimes can be found in the incoherent dynamics: A 
regime where the mean vibrational excitation is increased in each cycle, and a trap- 
ping regime for populations initially below a certain vibrational quantum number 
n where the condition w- mc (n + l\n ) = holds. The latter we will call, in the 
following, incoherent trapping states since they emerge from a dynamics bare of 
any coherence. 

In fact, incoherent trapping states emerge from a breakdown of the ion-laser cou- 
pling and are thus related to the second type of trapping states with coherent dy- 
namics discussed before. To see this relation, let us consider the condition for an 
incoherent trapping state at vibrational quantum number n : 

w iac (n +l\n ) = = \ {1 - exp [-7n o 0?) r s ]} , (45) 

which is equivalent to the condition 

7n (") = 0. (46) 
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From their definition (33) follows 7n ( ? ?) °c [/(^o! 7 ?)] 2 ? such that the condition for 
incoherent trapping states (46) is equivalent to the condition of nonlinear coherent 
trapping states (21). This is because both types of trapping states rely on the same 
mechanism, namely, the breakdown of the laser-ion coupling strength due to the 
vanishing of Franck-Condon type overlap integrals of vibrational wave functions. 
Thus, we deal with a trapping mechanism that appears both in the coherent and the 
incoherent regime and therefore does not rely on coherent dynamics. 

Incoherent trapping states may be conveniently prepared with no need for coher- 
ence. Decoherence effects will not prevent the emergence of these trapping states. 
For small vibrational quantum numbers one typically would need only few cycles 
such that little decoherence is accumulated. However, for generating highly excited 
vibrational number states, the mechanism of incoherent trapping states is superior. 
Then a rather large number of cycles is required, where in each cycle decoherence 
effects would eventually prevent any coherent mechanism from working. 

Let us now study the dynamic increase of the vibrational excitation in the incoher- 
ent regime. For large saturation s and/or large interaction times t s , the transition 
rates w- mc can be approximated as 



iyino(ra+l|ra) « - 



1 ^n,no(ii) 
no (v) 



w- mc (n\n) « - 1 + ^ £n,no(»?) 



(47) 
(48) 



where 5 n m denotes the Kronecker delta symbol and the sum extends over all pos- 
sible trapping numbers n (r)) at the given Lamb-Dicke parameter 77. This form 
of the transition rates is an excellent approximation for a wide range of parame- 
ters and holds for the following examples. Thus, in regions of vibrational numbers 
n=^n (r]), the transition rates simplify to 



[n + l\n) 



w iQC {n\n) « -, 



(49) 



and the recurrence relation (42) for the incoherent dynamics reduces to 

Pn(tk+l) = \ P n (tk) + \ Pn-l(tk)- 



(50) 



Due to the constant transition rates, now the evolution towards a trapping state is 
not slowed down by quasi-trapping conditions as may be found in the coherent 
rates. 

Using Eq. (50) it can be shown that, far from any trapping state, the mean vibra- 
tional excitation and variance obey the following dynamics: 
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(n(t k+1 )) = (h(t k )) + i, 
([An(t fc+1 )f> = ([An(t k )} 2 ) + \. 



(51) 
(52) 



These relations show that in each cycle, on the average, half a vibrational quantum 
is created in the ionic motion. Moreover, the mean excitation increases faster than 
the variance so that for a large number of cycles k the relative variance ( An 2 ) / (n) 
reaches the stationary value, 



A statistics of this type corresponds to a noise level half the shot-noise limit and is 
a signature for amplitude squeezing. In fact, the solution of Eq. (50) can be given 
as a sum of Binomial distributions, each with vibrational mean number l+k/2 and 
variance k/A, that are weighted by the initial number statistics Pi(t ): 



Near the trapping state \n ), however, the relative variance deviates from the sta- 
tionary value (53), decreases below 1/2 and reaches zero when the trapping state is 
finally populated with unit probability. 

Clearly, this behaviour is idealised by the fact that spontaneous recoil effects during 
the periods of optical pumping have been neglected in the analysis. However, using 
quantum-trajectory methods we have included these effects consistently to obtain 
numerical results. These results, that will be shown in the following, are based on 
the parameters ji/u — 9.5, 72/^ = 3.3, 771 ~ 772 = 0.142, and |f2 p | 2 /7 2 = 5.0 for the 
optical pumping processes. 

In Fig. 3 the numerical results for the relative variance of vibrational quanta is 
shown in dependence on the number of cycles for the incoherent trapping state 
I no = 50) for varying effective saturations 7sr s . It can be observed that for the 
larger values of 7st s the relative variance transiently tends to the value 0.5, but 
then increases due to the photon scattering until the trapping state is approached. 
There the relative variance decreases even below 0.5, but only as a local minimum, 
to monotonically increase afterwards. Furthermore, for low "saturation" (see dotted 
curve) the Binomial regime is not even reached as an initial transient behaviour. 

These results show that due to the recoil effects in the optical pumping the trapping 
states are not truly stationary. Instead the vibrational populations cross the trapping 
state to later be partly accumulated at approximately equidistant further trapping 
states at n ps 100, 150, 200, . . .. However, it should be kept in mind that the rela- 
tive variance is very sensitive to minimal populations at high vibrational quantum 
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Fig. 3. Relative variance versus the number of cycles k for the incoherent trapping state 
| no = 50) with r] = 0.268. Spontaneous recoil effects during the optical pumping are in- 
cluded; saturation parameter is 7sr s = 0.2 (dotted), 1.0 (solid), 2.0 (dashed). 

numbers. Thus a low percentage of population distributed around or higher than 
the trapping state drastically increases the relative variance. In this sense a more 
suitable property is the maximum population reached in the trapping state. 

An example is shown in Fig. 4 where vibrational statistics is shown at different 
stages of the evolution for the same trapping state |50). In fact it can be seen that 
the higher the value of 7sr s the higher the population in the trapping state at cy- 
cle k = 200. Clearly in the further time evolution this population decreases again 
by scattering events that eventually distribute the population at higher vibrational 
levels. Nevertheless, if the laser interactions are stopped appropriately when a high 
population in the trapping state is reached, maximum populations up to 60% at the 
level n = 50 can be obtained. 

It is noteworthy that for decreasing Lamb-Dicke parameters the trapping states 
have increasing energies. In Fig. 5 we show certain combinations of Lamb-Dicke 
parameters and vibrational quantum numbers n , that satisfy trapping. These are 
directly obtained from condition (46), together with Eq. (34). 

As to practical implementation of the incoherent method, experiments on a single 
ion have been performed that indicate breakdown of the ion-laser coupling at small 
values of the Lamb-Dicke parameter and rather highly excited trap levels [6,20]. In 
these experiments certain metastable states of the ion's vibrational amplitude with 
5-30/im spatial extension within the 400yum-wide focus of saturating laser light 
have been observed. In steady state, one of these metastable states was found oc- 
cupied by the unperturbed, moderately laser-cooled ion. The ion jumped to another 
one upon external perturbation: small variation of the cooling power by minute 
laser detuning, irradiation by a weak pulse of resonant radio frequency (l-2MHz, 
typically), etc. Each stable amplitude corresponds to a number of mean vibrational 
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Fig. 4. Incoherent dynamics, including spontaneous recoil effects, of vibrational pop- 
ulations for k = 0, 50, 100, 150, 200 cycles (gray scales from black to white) and 
Lamb-Dicke parameter rj = 0.268 with trapping state \uq = 50). Other parameters are 
jst s = 1(a), 10(b), 1000 (c). 

excitation that is associated with a particular level, a trapping level, where the laser- 
ion coupling breaks down for the given small Lamb-Dicke parameter. Thus the 
nonlinearity in the laser-ion coupling as discussed here, is of great importance even 
with a small Lamb-Dicke parameter r] that makes trapping states appear at very 
large vibrational quantum numbers. These observations demonstrate that an indi- 
vidual ion confined in a harmonic potential well indeed may be prepared in a highly 
excited incoherent trapping state. 



6 Summary and conclusions 



We have studied ways of generating various kinds of "trapping states" of the quan- 
tised centre-of-mass motion of a trapped ion. The first kind is of the micro-maser 
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Fig. 5. Pairs of values for the trapping state number no and the Lamb-Dicke parameter n 
that obey the nonlinear trapping-state condition [Eq. (21)]. 

type. In this case, the ion, at a particular vibrational quantum number, undergoes 
a complete Rabi cycle of excitation, such that higher quantum states cannot be- 
come excited. This mechanism, however, is highly fragile with respect to any kind 
of decoherence. Moreover, when one intends to prepare a certain vibrational num- 
ber state, quasi-trapping conditions may exist for other (lower) number states. This 
feature leads to a substantial increase of the time needed for the preparation. 

Another way of implementing a trapping-state condition is based on the nonlinear- 
ities that appear in the dependence of the atom-field coupling on the vibrational 
excitation. These nonlinearities arise as Franck-Condon-type overlap integrals due 
to the exchange of momentum between the ion's centre-of-mass motion and the 
absorbed and emitted laser photons. At certain vibrational quantum numbers the 
nonlinearly modified interaction strength vanishes, leading to the expected trap- 
ping states in a way quite different from the trapping mechanism based on the 
completion of Rabi cycles. This trapping mechanism, when used in its coherent 
version, reduces the preparation time, although this method is still highly sensitive 
with respect to decoherence. 

The particular advantage of the nonlinearity-based trapping condition consists in 
the vanishing coupling for certain vibrational excitations. This feature can be ex- 
ploited for the generation of trapping states even in situations where substantial 
decoherence prevails, since now the trapping state is solely determined by the se- 
lected Lamb-Dicke parameter. Moreover, the transition rates for approaching the 
trapping state are more or less independent of the vibrational excitation, which 
leads to substantial decrease of the time needed to prepare a desired trapping state. 
Thus, one avoids the sensitivity to decoherence, especially when preparing vibra- 
tional number states of large quantum numbers. 
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